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Abstract 

The realistic equation of state of strongly interacting matter, that has been success- 
fully applied in the recent hydrodynamic studies of hadron production in relativistic 
heavy- ion collisions at RHIC, is used in the Friedmann equation to determine the 
precise time evolution of thermodynamic parameters in the early Universe. A com- 
parison with the results obtained with simple ideal-gas equations of state is made. 
The realistic equation of state describes a crossover rather than the first-order phase 
transition between the quark-gluon plasma and hadronic matter. Our numerical cal- 
culations show that small inhomogeneities of strongly interacting matter in the early 
Universe are moderately damped during such crossover. 
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1 Introduction 

One expects that the studies of relativistic heavy-ion collisions may bring us 
information about the properties of matter in the early Universe (at the time 
when its age was about 10 microseconds and its temperature was about 200 
MeV). In this paper we give an explicit example showing how such expec- 
tations may be realized in practice. Our main point is that the theoretical 
analysis of the hadronic observables collected at RHIC (Relativistic Heavy 
Ion Collider at BNL) pinned down the form of the equation of state (EOS) 
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of strongly interacting matter. The use of this EOS in the Friedmann equa- 
tion allows us to determine the precise time evolution of the thermodynamic 
parameters in the early Universe. 

Until recently, the soft hadronic observables measured at RHIC have been 
described with the help of the hydrodynamic models which incorporated the 
first-order phase transition between the quark-gluon plasma and hadronic mat- 
ter [T|2] . The main success of this approach was the correct reproduction of 
the hadron transverse- momentum spectra and the elliptic- flow coefficient t> 2 . 
On the other hand, the hydrodynamic models using the first-order phase tran- 
sition had severe problems with the reproduction of pionic correlation func- 
tions. The latter contain the information about the space-time dimensions of 
hadronic systems formed in the collisions. Thus, the evident successes in the 
reproduction of the hadron momentum distributions were obscured by the fail- 
ures connected with the analysis of the space-time dimensions of the produced 
system. 

The large part of those difficulties, known in the literature as the HBT puz- 
zle, was removed by the use of the EOS with the crossover phase transition 
01]. The crossover phase transition is not a genuine phase transition — all 
thermodynamic parameters change smoothly as functions of the temperature, 
however, such changes may be quite sudden in a narrow temperature range. 
The crossover phase transition without a distinct soft point (i.e., without a 
thermodynamic region where the ratio of pressure and energy density is very 
small and the sound velocity possibly drops to zero) leads to a faster time 
development of the system and helps to reproduce the correlation data. 

The use of the simple bag EOS (leading to the first-order phase transition) in 
the hydrodynamic equations describing heavy-ion collisions is very much rem- 
iniscent of the situation in the cosmological calculations where, in most cases, 
the strongly-interacting matter is treated in the same naive way. Therefore, 
it is interesting to search for the effects of using the realistic QCD EOS in 
the Friedmann equation. In particular, the results of the lattice simulations of 
QCD with physical quark masses, which suggest the crossover rather than the 
first-order phase transition, should be incorporated in such analyses. The work 
in this direction has been recently done in Ref. [5] where the lattice results 
from [6J were used. 

Our present approach also uses the lattice results but differs from Ref. [5] in 
several important aspects. Firstly, we use different lattice results as an input 
[7]. Secondly, our EOS is not based exclusively on the lattice results. At lower 
temperatures we use the realistic hadron-gas model, which turned out to be 
very successful in the thermodynamic analyses of hadron production in heavy- 
ion collisions |8J. The results of the hadron-gas model at lower temperatures 
are combined with the lattice results at higher temperatures in the thermody- 
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namically consistent way described in [9] . However, the most attractive feature 
of our EOS is the fact that it was successfully implemented in 2+1 [TO] and 
3+1 [11] hydrodynamic codes describing heavy-ion collisions. Thus, in our case 
the knowledge from heavy-ion physics is transferred directly to cosmology. We 
note that an EOS sharing very similar features to ours has been very recently 
constructed in Ref. [12] . 

In the second part of the paper we consider the time development of the 
energy density fluctuations using the theory developed by Schmid, Schwarz, 
and Widerin [T3JT3]. By doing explicit numerical calculations we show that 
small inhomogeneities of the strongly interacting matter are moderately (by 
about 30%) damped during the realistic crossover transition. 

Altogether, our results indicate that there are small chances for observation 
exotic phenomena connected usually with the first order phase transitions (the 
latter are discussed, for example, in [T5]). 



2 Equations of state for strongly interacting matter 

The QCD EOS used in this paper was introduced and discussed in greater 
detail in Ref. [9]. As mentioned above, it interpolates between the hadron- 
gas results and the lattice results. All thermodynamic properties follow from 
the temperature dependence of the sound velocity c 2 s depicted in Fig. [1] The 
smooth dependence of the sound velocity on the temperature indicates that 
there is no genuine phase transition in the system. On the other hand, c 2 (T) 
has a characteristic minimum at T ~ 170 MeV. This behavior reflects the 
presence of the sharp crossover in this region, i.e., the region where the energy 
density increases strongly within a narrow range of temperature. We stress 
that the values of c 2 . differ from 1/3. This indicates that the matter under 
consideration does not behave like an ideal gas of massless particles. 

Various thermodynamic quantities characterizing our QCD EOS may be ob- 
tained from the thermodynamic identities valid at zero baryon chemical po- 
tential 

e + P = Ta, de = Tda, dP = adT, c 2 = ^ = (1) 

de Tda 

Here s, P, T, and a are the energy density, pressure, temperature, and entropy 
density, respectively. In Fig. [2] we show the energy density scaled by T 4 (solid 
line), that was obtained directly from the function c 2 (T) depicted in Fig. [TJ 

For comparison, we also show the energy density Em used by the Milan group 
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Fig. 1. Solid line: temperature dependence of the sound velocity for the realistic 
QCD EOS constructed in [9] and used in this paper. Dashed line: the limit for 
massless and non-interacting particles. 

in Ref. [5J. The function e M (T) is again scaled by T 4 (dashed line). We ob- 
serve that our values of e/T 4 are larger at lower temperatures and smaller 
at higher temperatures. Larger values of e/T 4 at lower temperatures are con- 
nected with the inclusion of all known resonances in the hadron-gas model. 
On the other hand, smaller values of e/T 4 at higher temperatures indicate 
stronger interactions in the plasma phase. 

In order to compare the results obtained with the realistic QCD EOS with 
the results obtained with the bag EOS we introduce the following definitions 




71-2 4 ti-2 

e qgp = 9qgp~^T + B, -Pqgp = g q gp—T — B, if T > T c , (2) 



7T 



e =a —T 4 
* ^30 ' 



7T 



if T <T r 



(3) 



Here B is the bag constant and T c is the critical temperature defined by the 
formula 



905 



1/4 



.(0qgp -9*W 



(4) 



Equation (|2|) describes the energy density and pressure of the weakly inter- 
acting quark-gluon plasma, whereas Eq. (EJ) describes the energy density and 
pressure of the massless pion gas. For two quark flavors g qgp = 37, and for 
three flavors g qgp = 47.5. 

In the numerical calculations presented below we use g qgp = 37 and consider 
two values of the bag constant: B = (235 MeV) 4 and B = (165 Me V) 4 . 
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Fig. 2. Solid thick line: the energy density e obtained from the function c 2 (T) shown 
in Fig. [Hand scaled by T 4 . Dashed thick line: the energy density em used in Ref. 
[5], also scaled by T 4 . Thin dashed lines: the Stefan-Boltzmann limits for the ratio 
e/T 4 for the bag EOS with 2 (lower line) and 3 (upper line) flavors. 

They correspond to two values of the critical temperature: T c = 169 MeV 
and T c = 119 MeV. 

We stress that the use of the bag EOS as a reference point for the calculations 
with the realistic QCD EOS is ambiguous. For example, we easily find that the 
bag EOS yields the Stefan-Boltzmann limit £ qgp /T 4 = 377r 2 /30 for two flavors 
and £ qgp /T 4 = 47.57r 2 /30 for three flavors. Figure [2] shows that the realistic 
value of e/T A does not reach any of those limits. The difference indicates that 
the interactions in the plasma are non-negligible at the temperatures as high 
as 1 GeV. 

Below, we consider the time evolution of the thermodynamic quantities gov- 
erned by the Friedmann equation. The initial time corresponds to the moment 
when T = T = 500 MeV and e(T ) = 100 GeV/fm 3 . This value is close to the 
energy density obtained at the same temperature with the bag EOS with two 
flavors. This coincidence suggests that it is useful to consider the bag EOS 
with two flavors as a reference EOS for T < 500 MeV. 



3 Inclusion of the electro-weak sector 



In the temperature range considered by us in the Friedmann equation, one 
should include also the energy density and pressure of the particles from the 
electro-weak sector. Typically they are treated as massless particles, whose 
effective number of degrees of freedom is g cw = 14.45 [16J. Thus we write 
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Fig. 3. Temperature dependence of the total energy density er = e + e ew and the 
reference energy density er obtained with the bag EOS. 



£R — £ + £ew> Pr — P + -Pew, (5) 

where Er and Pr are the total energy density and pressure used in the Fried- 
mann equation and 



fi'cw gg J- > -*ew fi'cw gg J- ? C"cw ^rp ■ [y) 



In the similar way we construct the modified bag EOS 



2 2 

e B = g ^-T 4 + B 1 P B = gi ^T 4 -B, if T>T C , 

e B = g 2 K \\ p B = g2 K l T \ if T<T C , (7) 



where g\ = 37 + g ow = 51.25 and #2 = 3 + g cw = 17.25. 

Since g qgp — g n — g\ — g 2 , the value of the critical temperature for fixed B re- 
mains the same as in the case without the electro-weak sector. We note that we 
use now the same subscript to denote the thermodynamic variables above and 
below T c . The temperature dependence of the functions £r(T) and e B {T) in 
the temperature range 0.2 GeV <T< 0.5 GeV is shown in Fig.|3j We observe 
that the two functions converge in the region 0.45 GeV < T < 0.50 GeV. 
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4 Solutions of the Friedman equation 



Assuming that the expansion of the Universe is isentropic, the Friedmann 
equation in our case takes the form 



de 



R 



dt 



In the next step we use Eq. (JHJ) and switch to the temperature as the inde- 
pendent variable. In this way we get 



3cr 



dT, 



R 



l8irG(e + e e 



dt 



(e + e ew + P + P ew ). 



(9) 



We note that c~ 2 , a, e, and P in Eq. (Q characterize the realistic QCD EOS, 
while cx ew , e ew , and P ew are defined by Eq. fl§]). The temperature obtained 
from Eq. has a subscript R in order to distinguish it from the temperature 
profile T B obtained for the bag EOS 0. 

The initial condition for solving Eq. ([9]) has the form 

T R (t ) = 500 MeV. (10) 

The corresponding initial energy density 

e R ,o = e(T R (t )) + e cw (T R (t )) = 138 GeV/fm 3 . (11) 



Assuming that the evolution for t < to is dominated by matter with the 
radiation-like EOS we set 

^ (i2) 



The numerical solution of Eq. (J5J), i.e., the function T R {t) as well as the func- 
tions e R (t) and Pn{t) are shown in Figs. H] and [5] (thick solid lines). The results 
obtained with the realistic QCD EOS are compared with the results obtained 
with the bag EOS (thick dashed lines) and the same initial conditions. The two 
figures differ by the value of the bag constant used in the reference calculation: 
B = (235 MeV) 4 and B = (165 MeV) 4 in Figs. 1 and EJ respectively. 

1 The reader may connect the subscript R with the realistic or RHIC equation of 
state. 
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Fig. 4. The energy density er (a), temperature Tr (b), and pressure Pr (c), ob- 
tained by solving the Friedmann equation with the realistic QCD EOS with the 
initial condition Tr = 500 MeV at t$ = 1.35 fis (solid lines). The dashed lines show 
the results obtained for the bag EOS with B= (235 MeV) 4 and the same initial 
condition, Tr^q) = Tr^q). The initial and final times for the reference first-order 
phase transition are denoted by t\ and t2, respectively. 

The time changes of the functions £_r(£) and £#(£) shown in Fig. H] are very 
much similar. On the other hand, the changes of the temperature and pressure 
are quite different in the region where the system governed by the bag EOS 
passes through the phase transition. For large values of time, the functions Tr 
and Pr become again close to Tr and Pr- 
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Fig. 5. The same as Fig. U but with the reference bag calculation done for 
B= (165 MeV) 4 . 



The value of the bag constant used in the calculation shown in Fig. H] implies 
T c = 169 MeV. This value agrees with the critical temperature used in the 
realistic QCD EOS. Strictly speaking, the realistic crossover phase transition 
does not allow for the unique definition of the critical temperature. Neverthe- 
less, in our approach we associate T c of the crossover with the local minimum 
of the sound velocity function shown in Fig. [TJ which is situated close to 
T = 170 MeV. 
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Fig. 6. Time evolution of the scale factor a(t) (solid line) compared with two pow- 
er-law expressions: (t/io) ' 50 (lower dashed line) and (i/io) ' 52 (upper dashed line). 

In Fig. Owe compare the same results obtained with the realistic QCD EOS 
with the reference calculations based on the bag model with B = (165 MeV) 4 . 
In this case the value of the bag constant corresponds to T c = 119 MeV. 
Clearly, the agreement between the two results is much poorer. 



5 Scale factor and conformal time 

The time evolution of the scale factor a(t) is given by the formula 



In Fig. [6] we show the time dependence of a(t)/a(to) obtained with our energy 
density function 6r (solid line). For comparison, we also show the two power- 
law functions (t/to) p , where p = 0.50 (lower dashed line) and p = 0.52 (upper 
dashed line). Clearly, the use of the realistic EOS changes slightly the standard 
massless-gas behavior corresponding to p = 0.5. 

In the following we use the approximation of a(t)/a(to) by the function of the 
form (t/t ) p and in the numerical calculations we set p = 0.5170. In this case, 



2 The value p = 0.517 follows from the power-law fit done in the whole time in- 
terval between t = 1.35 /us and t = 220 /us. For t < 10 fis the function a{t)/a{to) is 
much better described by the formula (t/to) ' 5 , however, the presence of the phase 
transition modifies the effective power-law behavior at larger times. We note that 
the change of p from 0.5 to 0.517 results in the 12% change of a at t = 220 /xs. 




(13) 
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the Hubble constant is 



"=4=f < 14 > 

a at t 



For further applications it is useful to introduce the conformal time rj defined 
by the relation 



We further define the parameter k and the variable x 
- , x = k r} = - 



(15) 



dr] 1 
dt a' 

and the quantity H that is related to the Hubble constant, 

U = aH = ^-. (16) 
a drj 

Simple calculations lead to 

(l-p)a {l-p)r] 



a = , x = k rj = — . (1J 



6 Cosmological inhomogeneities at the QCD phase transition 



The character of the QCD phase transition affects not only the time depen- 
dence of average thermodynamic quantities but it has also influence on the 
primordial density fluctuations in the early Universe. If the QCD phase tran- 
sition was of the first order, with the sound velocity vanishing at the criti- 
cal temperature, the density fluctuations would develop large peaks and dips 
[T3fTi] . For the crossover transition we do not expect any spectacular phenom- 
ena, however, it is interesting to analyze the development of fluctuations in 
the realistic background analyzed in the previous Sections and to check if the 
drop of the sound velocity at T ~ 170 MeV may have any noticeable effects 
on the inhomogeneities. 

Below we study the system of equations introduced by Schmid, Schwarz, and 
Widerin [I5II3] . They read 
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±6' + 3(c 2 s - w)5 = - 3(1 + w) a, (19) 

^C=^Vw-4a, (20) 

Ltf + (l _ 3 ^ = -c 2 A<5 - (1 + «,) A a , (21) 
1 k 4k 



a=-|(l + 3<4) fe. (23) 



Equations ffl~9]) and (1201) follow from the energy-momentum conservation, Eqs. 
fl2TT) and (122]) from the 3-divergence of the Euler equation of general relativity, 
and Eq. (1231) from the Einstein i?Q-equation. The prime denotes the derivative 
with respect to the conformal time rj. The quantities 5 describe the energy 
density fluctuations (density contrasts), ip's are related to the fluid velocities 
(peculiar velocities), and a defines the correction to the temporal part of the 
metric tensor (lapse function). The function w is the ratio of pressure and 
energy density for the background evolution, and c s is the sound velocity of 
the background (in our case we use here the results of solving the Friedmann 
equation presented in the previous Sections). 

Equations flU and fl5T| refer to the QCD matter, whereas Eqs. ([20]) and ([22]) 
refer to the electro-weak sector. On the other hand, Eq. ([2"B"]) refers to the two 
components, hence, in this case we have 



w R 



^sR 



Pr_ 
dP R 



ds 



R 



P + P 

1 ~ 1 ew 

£ 1 ^ew 

dP + dP e , 
de + cfeew 



a + a c 



c s 2 a + 30c 



(24) 



and 



5e R 5e + 5e ew 5 + 5 ew (e cw /e) 
°R = = ; = i,/ 7"T~- ( 25 ) 

£R £ + £ew 1 + l^cw/^J 

The functions P, P ew , e, e ew , a and a ew are known if the temperature depen- 
dence on time is established by Eq. (J9j) . 

The parameter k in Eqs. (jT9]) - (l2"3]) is the comoving wavelength, k = & p h ys a - in 
the numerical calculations we switch from the conformal time 77 to the variable 
x = koT). In this case, using the variable £ = k/ko, see Eq. ( 1T8]) . we may write 
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1 ■ 1 — p dd 1 — p dS ■ 

—5 = 7]— = x— = qxd, 

ti p dr/ p dx 

— = - — -r]k = q^x. (26) 
ti p 

The dot in f )26|) denotes differentiation with respect to x. Using substitutions 
of the type f )26|) in Eqs. f|T9|) -f f23|) we obtain a system of 5 coupled ordinary 
differential equations for 5 unknown functions: 5, 5 ew , if), ip ew , and a: 



qx S + 3(Cg — w) 5 = q£xip — 3(1 + w) a, 
qx 5 ew = q£x ip ew -4 a, 



qx ip + (1 — 3w) if) = —q^x 



c 2 J+ (1 +w)a 



Vw = -g (4w + 4a) 



(q£ x ) 2 + ^(l +WR ) 



a = -|(l + 3c2 H ) (5 K . (27) 



We analyze the evolution of inhomogeneities in the time interval between 
t = 1.35 /is and t = 220 /xs. Since x = (t/to)^ 1_p \ the starting value of x is 1 
and the final value of x is about 12. During this evolution, the temperature 
of the system drops to 53 MeV, thus we explore the whole temperature range 
where the sound velocity has a minimum, see Fig. [TJ 

In the numerical calculations we consider the subhorizon modes that satisfy 
the condition k/T-L ^> 10- This condition is equivalent to qxt; ^> 1 or simply 
to £ 3> 1. In Fig. [7] we show our results obtained for £ = 10 with the initial 
conditions: 5 = 1.0, if) = ip ew = 0, and 6 ew = —0.5 (upper part) and 5 CW = 
(lower part). The density contrast of the strongly interacting matter 5 (solid 
line) and of the electro-weak sector <5 ew (dashed line) oscillate with frequencies 
determined by the value of the sound velocity in medium. This is indicated 
by the thin solid line representing the function 



cos 



£ jcl(x')dx'\ . (28) 



The frequency of 5 is reasonably well reproduced by Eq. (|28|) . however, in the 
considered time interval the amplitude of 5 decreases by about 30%. 

In Ref. [H] a schematic crossover transition between the ideal quark-gluon 
plasma and a gas of massless hadrons was analyzed (in addition to the much 

3 The superhorizon modes, with k/H <C 1, are independent of the thermodynamic 
behavior. The density contrasts grow quadratically with k/H [13] . 
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Fig. 7. Time evolution of the density contrasts 5 (solid line) and 5 ew (dashed line) 
for £ = 10. The values of the initial conditions are given in the Figures. 

more elaborated case of the first order phase transition). A similar damping 
was observed in this case, that was identified with the reduction of the number 
of degrees of freedom during the crossover transition. During the first order 
phase transition this type of damping is also present but it merely reduces the 
huge amplification of the modes triggered by vanishing of the sound velocity. 

For the realistic QCD EOS the sound velocity does not drop to zero at T c and 
we deal only with the moderate damping of the inhomogeneities of strongly 
interacting matter. On the other hand, the inhomogeneities in the electro- 
weak sector, where the number of degrees of freedom is fixed, are not damped. 
They are well described by the function of the form fl28|) with a suitably fixed 
phase and c 2 s = 1/3 (not shown). 

The lower part of Fig. [7] differs from the upper part by the initial value of <5 ew 
that is set equal to zero in this case. Interestingly, the initial non-zero value of 
S triggers oscillation in the electro- weak sector. Similarly to the upper part, 
the amplitude of the oscillations of 5 is damped by about 30% during the 
considered time interval. 

In Fig. M we show our results obtained for £ = 5. Similarly to the previous 
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Fig. 8. The same as Fig. [7] but for £ = 5. 

case with £ = 10 we observe again a slow damping of the oscillations of the 
QCD matter. 

7 Conclusions 

In this paper we have considered the effect of the realistic QCD EOS on the 
evolution of the early Universe. We have compared the solutions of the Fried- 
mann equation with those obtained with the bag EOS. The use of the realistic 
QCD EOS has small effects on the time evolution of the energy density, how- 
ever, there are noticeable differences in the time evolution of pressure and 
temperature, which are more affected by the first order phase transition. Such 
difference may have an impact on various physical quantities if the cosmolog- 
ical measurements are done with high precision. The differences between the 
realistic QCD EOS and the bag EOS are the smallest if the critical temper- 
ature of the bag model is fitted to reproduce the critical temperature of the 
realistic crossover transition. 

In the second part of the paper we have considered the time development of 
the energy density fluctuations. We have shown that small inhomogeneities 




f = 5, <5(x )=l, (5 ew (xo)=-0.5, (/r(x )=l/'ew(xo)=0 




f = 5, 6(x )=l, <5 ew (*o)= tA(jr )=i/'ewUo)=0 
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of the strongly interacting matter are moderately (by about 30%) damped 
during the realistic crossover transition. Our results indicate that there are 
small chances for observation exotic phenomena connected with the first order 
phase transitions. 

This research was supported in part by the Polish Ministry of Science and 
Higher Education, grant No. N N202 263438. 
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